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TECHNICAL NOTE 4320 

COMPRESSIBLE LAMINAR PLOW AND HEAT TRANSFER ABOOT 
A ROTATING ISOTHERMAL DISK 
By Simon Ostrach and Philip R. Thornton 


SUMMARY 

The flow and heat transfer about a rotating isothermal disk are re- 
examined to include the effects of compressibility and property varia- 
tions. If viscous dissipation is neglected,, the compressible problem is 
correlated to the incompressible problem by assuming linear variations 
of viscosity and thermal conductivity with temperature. Certain in- 
accuracies in several previous incompressible solutions are noted and 
corrected herein. The effect of compressibility appears as a distortion 
of the normal coordinate and normal velocity component and as a multi- 
plicative factor in the heat- transfer coefficient, the Nusselt number, 
and in the expressions for the skin- friction components and torque re- 
quired to rotate the disk. 


INTRODUCTION 

The steady laminar motion of an incompressible viscous fluid about 
an infinite rotating disk was considered by von Karmen (ref. l) . The 
Navier-Stokes equations were reduced to ordinary differential equations 
by separation of variables, and these were then solved by the K&rm&t- 
Pohlhausen integral method. Cochran (ref. 2) corrected several errors 
that he noted in von Kfirm&i's solutions, and using the corrected inte- 
gral solutions as a first approximation he obtained more accurate re- 
sults by numerically integrating the differential equations. The flow 
was shown to be similar to that about a centrifugal fan: The fluid 

moves radially outward, especially near the disk, and to preserve con- 
tinuity there Is an axial flow toward the disk. 

The heat transfer from a uniformly heated rotating disk was first 
considered by Wagner (ref. 3) who used von K&rm&n 1 s uncorrected results 
to solve the energy equation neglecting dissipation and thereby derived 
a heat-transfer coefficient. In reference 4 the heat-transfer problem 
is also treated, but there Cochran's flow solutions are used and the 
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energy equation including viscous dissipation is solved. However , the 
analyses (refs. 3 and 4) are restricted to apply to very small rotation 
and very small temperature differences. For large rotations or heating 
one might anticipate that the compressibility and property variations of 
the fluid would be important. Therefore, the problem is reexamined here- 
in to include these effects. 


BASIC EQUATIONS AND BOUNDARY CONDITIONS ' 

The configuration to be studied is that of an infinite circular disk 
in the R,0 plane which is maintained at a uniform surface temperature 
and can rotate about the Z-axis (see fig. l). The equations in cylindri- 
cal coordinates expressing conservation of mass, momentum, and energy for 
an axially symmetric compressible viscous flow are, respectively, 


s m ^ pRU > + ^ “ 0 (la) 



where the symbols are all defined in appendix A and where 
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x> 


For a disk rotating with an angular velocity ft, the no- slip condi- 
tion for viscous fluids requires that 


U(R,0) = W(R,0) = Oj V(R,0) = ftR 
At large distances frcm the disk it is specified that 

U(R,«) = V(R,«0 = 0 

The thermal boundary conditions for an isothermal disk are 

T(R,0) = T w , T(R,“) = Too 


( 2 ) 

( 3 ) 

( 4 ) 


In order to nondimensionalize the basic equations and reduce than 
to ordinary differential equations, let 


R =» (b^/ft) 1 / 2 r, Z = (b^/ft) 1 / 2 z 

U = (ftbj 1 / 2 r f(z), V - (ftoj 1 / 2 r g(z), W - (ftoj 1 / 2 h(z) 
P = PoobooOrt ( z ) j p = Poo p; t(z) = (T - Tj/(T W - Tj 
and assume the viscosity and conductivity laws 

p = CWii^fp) k =« ^ ^ ^ C(z) 


( 5 ) 


Equations (l) then become 


^ (pr 2 f ) + ^ (prh) = 0 


p(f 2 + hf' - g 2 ) = 


a 

3TI 


Pr p hr 


tH 

p(2?g + hg') = ^ gj 

FhH' »-«'+! C h ’ - f )j + I r2f ’] 

^ 8 fc(z) .1 Pr pft 

- Cp(T w - Too) ^[f T j + ^ 


3 ( f “ h‘) 


( 6 ) 

(7a) 

(7b) 

(7c) 

(7d). 

.\2 


(-) 


where the primes denote differentiations with respect to z and where 
Cp/Pr is assumed to be constant. 
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The boundary conditions (eqs. (2) to (4)) become 


«*\ 

o 

11 

o 

& 

II 

o' 

•w. 

'tH 

g(0) = 1 

(8) 

f(°°) = g(“) = o 

(9) 

t(0) = 1; 

t(oo) sa 0 

(10) 


If equations (7a) to (7c) were written in terms of the velocities 
and their derivatives , they would be recognized as the three-dimensional 
boundary- layer equations for compressible axisymmetric flow. The boundary- 
layer nature of the flow is also pointed out in references 1 and 4. Equa- 
tion (7d) is essentially an equation for the normal pressure gradient that 
is negligible within the boundary- layer approximation. The pressure in 
this problem is therefore constant. 

In order that equations (7a) to (7c) and (7e) be a consistent set, 
it is necessary that the last term of (7e) be neglected because of its 
r- variation. To find the conditions under which this term can be omitted, 
it is noted that the coefficient appearing in equation (7e) Is essentially 
given by (y - l)Pr 0 jM 2 /Re{^(T w /T 00 ) - lj . Since the r is not of unit order 

but actually equal to the square root of the Reynolds number Re, the last 
term of equation (7e) can be neglected if 

(r - DPrJl 2 

7p « 1 ( 11 ) 


This condition is obviously satisfied for relatively moderate rotations 
and large temperature ratios. The second- last term of equation (7e) is 
even smaller than the last terms and therefore is also negligible if in- 
equality (ll) Is satisfied. These terms are neglected in the remainder 
of the report. 


Equations (7) can be reduced to the incompressible equations for 
this problem as given essentially in reference 4 by taking C(z) to be 
a constant, C (implying a linear viscosity- temperature law) and letting 



T = S, rt = H= constant J 


( 12 ) 
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Equations (7) then become 


2F + H 1 =0 

(l3a) 

F 2 + HF’ - G 2 = F" 

(13b) 

2FG + HG' = G” 

(13c) 

11 

o 

(13d) 

PrHS 1 = S" 

(l3e) 


where here the primes denote differentiation with respect to £ . The 
boundary conditions are 


F(0) = H(0) = 0; G(0) = 1 

F(oo) = G(co) = 0 
S(0) « 1; S(oo) = 0 


(14) 

(15) 

(16) 


Equations (13a) to (13c) with the boundary conditions (eqs. (14) and (l5)) 
were solved in references 1 and 2. Equation (l3e) was solved in refer- 
ence 4 with the boundary conditions (eq. (16)), where S = Qj_ and 
Pr = ycr. Thus, once the solutions of equations (l3a) to (l3e) are known, 
the solutions can be transformed back to the compressible solution by 
means of equations (5), (6), and (12). The effect of compressibility, 
therefore, is to distort the normal coordinate and normal velocity. 


SOLUTIONS 

After F was eliminated by means of equation (13a), the dependent 
variables G and H were determined by numerical integration of equa- 
tions (13b) and (13c) along with the boundary conditions (l4) and (15) 
on an IBM 650 high-speed computing machine. The problem was treated as 
an initial- value problem, and a forward integration technique similar to 
that described in appendix A of reference 5 was employed. Five-point 
integration formulas were used throughout. The recalculation of these 
functions was made here to extend them to large values of the independ- 
ent variable. This extension is necessary for an accurate calculation 
of the temperature distribution and hence the heat transfer. 

The dependent variable S was then calculated ■ from equation (l3e) 
on a desk computer for Pr = 0.72 using the results for H as calcu- 
lated herein. 
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RESULTS 

Velocities 

From equations (5) and (12) it can be seen that the velocities are 
given essentially by the functions F, G, and H. The functions H, H', 
H" , G , and G' are listed in table I. H and G differ slightly from 
Cochran’s results (ref. 2) with H(<») tending to the value -0.8845 as 
compared to -0.886 for Cochran (ref. 2). The velocities, which are 
plotted in figure 2 against £, tend to their limits very rapidly, as 
they should for a boundary- layer flow. 


Temperatures 

From equations (5) and (12) it can be seen that the temperature is 
essentially equal to S. The function S is plotted in figure 3 against 
£ along with the function for a = 0.514 (Pr =* 0.72 for air) from 

reference 4. The function S, as calculated in the present report, is 
more accurate because the H function, is known to four decimal places 
for all £, while in reference 4 the values of H were taken from ref- 
erence 2, where H was tabulated for values of £ to 4.4. 


Heat Transfer 


The average heat-transfer coefficient is defined as 


h = 


/ ^ (si)z=c 


(IV - TJA 


>z=o 

A(T W - TJ 


dA 


However, because of equation (6), 


K p* As\ AA 1 / 2 
Vc p oo W$. 0 W) 


h = 


- Vc Ko 




(17) 


The constant C (frcm the viscosity and conductivity laws) can be chosen 
for best correlation. 


To compare with the previous work in which only small 
differences were considered, C can be taken to be unity, 
the present case 

, „ A/ 2 

h = 0.329 k*, * 




temperature 
Therefore, in 


( 18 ) 
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Wagner in reference 3, defining h as in equation (18), obtained 

L \ l/2 


h = 0.339 k 




(19) 


By means of an enthalpy balance, Wagner also found 


/«W ! 

b = 0.335 It. 


( 20 ) 


where in equations (l9) and (20) he assumed Pr = 0.74. The difference 
between the present results and Wagner’s (ref. 3) may be due to two 
things: (l) Wagner used von K^rmdn’s (ref. l) uncorrected velocity pro- 

files in his calculations, and (2) Wagner used a larger value of Pr 
than is used herein; this would tend to give a larger initial slope for 
the temperature function according to figure 7 of reference 4. If the 
results of reference 4 were used to determine h, the coefficient would 
be 0.286. The difference occurs because an incorrect energy equation 
was used in reference 4, which led the- parameter Pr/y (taken as 0.514 
there) to appear rather than the Prandtl number itself (see eq. (l3e)). 
This point was first noted in reference 6 and is discussed in detail in 
appendix B. 


The heat-transfer data of references 6 and 7, along with the analyses 
of references 3 and 4 and the present report, can al 1 be expressed in 
terms of an average Nusselt number defined as 



00 


(21a) 


where 



(21b) 


for the present report. In figure 4, the ratio of the Nusselt number to 
the square root of the Reynolds number from various analyses and experi- 
mental studies is plotted against the angular velocity, where, for the 
present analysis, it is assumed that C = 1. If the mean of the heat- 
transfer data as correlated in reference 6 is used, only the experimental 
results of reference 7 vary with ft. The data in reference 7 give re- 
sults that are much higher than any others, while the mean of the data 
for reference 6 are in good agreement with the present analysis. Refer- 
ence 6 points out that the results of reference 7 are too high because 
of heat losses through the insulation and that the results are also 
affected at low rotational speeds by free convection (due to gravitational 
force ) . The latter effect is clearly shown in figure 4 . The value of 



8 


' NACA TN 4320 


h, as given by equation (20), is used in equation (21a) to plot the anal- 
ysis of reference 3 in figure 4, while the value of h corresponding to 
Pr/Y = 0.514 is used for reference 4. The latter is used in order to 
illustrate the importance of using the correct form of the energy equa- 
tion. If the results for Pr/r = 0.72 from reference 4 are used, the 
value of Nu/-y/Re is 0.35. 


Skin Friction and Torque 

Because of angular symmetry and the form of the normal velocity 
component (eq. (5)), the expressions for the components of tangential and 
radial skin friction are 


Radial skin friction = J' ( T ZR ) Zf -o ^ = f Mv [fz)z=0 ^ ( 22a ) 

Tangential skin friction = / < 22t > 

where A equals the area of a finite disk of radius Rq. Thus, if equa- 
tions (5), (6), and (12) are used, these expressions can be written as 


Radial skin friction = | xpJl-JCoji R§F'(0) = 1.069 p Jl-y/GoJl R§ 

(23a) 

Tangential skin friction = R^G'(o) = -1.289 p^-/Cu^0 Rq 

(23b) 

The torque, or the rotational moment necessary to turn the disk, 
is 


Torque = f EdA . 2 n f ° (t^)^ K 2 ® (24a) 

Again, if equations (5), (6), and (12) are used, this equation becomes 
Torque = § pJl-^Cvji rJg’(O) = -0.98 p Jl-y/Ooji Rq (24b) 

Von Karmen (ref. l) derived a similar expression using the angular mo- 
mentum leaving the cylindrical surface formed by the disk and the gaseous 
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boundary layer. The coefficient in reference 1 was 0.92. Again, the 
results are affected by the compressibility as evidenced by the -Jq 
which appears in the equations. 


SUMMARY OP RESULTS 

The flow and heat transfer about a rotating isothermal disk have 
been reexamined to include the effects of compressibility and property 
variations. Por flows in which a relatively highly heated or cooled 
disk is rotating with a moderate velocity so that viscous dissipation is 
negligible, the compressible problem is correlated to the incompressible 
problem by assuming linear variations of viscosity and thermal conduc- 
tivity with temperature. Certain inaccuracies in several previous in- 
compressible solutions have been noted and corrected herein. The effect 
of compressibility appears as a distortion of the normal coordinate and 
normal velocity component and as a multiplicative factor in the heat- 
transfer coefficient, the Nusselt number, and in the expressions for the 
skin-friction components and torque required to rotate the disk. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 2, 1958 
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APPENDIX A 


SYMBOLS 

A 

surface area of disk 

C(z) 

function in viscosity and conductivity laws, eq.. (6) 

c p 

specific heat at constant pressure 

C v 

specific heat at constant volume 

r(g),G(6)#H(5) 

"incompressible" dimensionless velocity functions 

f(z),g(z),h(z) 

compressible dimensionless velocity functions 

E 

average heat-transfer coefficient, q/A(T w - T^) 

k 

coefficient of thermal conductivity 

M 

Mach number, 

Nu 

Nusselt number, hf^/k^ 

P 

pressure 

Pr 

Prandtl number, CpM./k 

%. 

dimensionless temperature function defined in ref. 4 
equal to S(£) 

<1 

total heat flux 

R,0 ,Z 

cylindrical coordinates, Z normal to disk surface 

R 

gas constant, P =» pRT 

Re 

Reynolds number, £2R 

T,Z 

dimensionless cylindrical c.oordinates 

s(0 

"incompressible" dimensionless temperature function 

T 

temperature 

u,v,w 

radial, tangential, and normal velocities 
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X 

5 

n 

o 

«,n 

p 

a 

r 

t zr 

T Z0 

$ 

ft 

Subscripts : 
w 


ratio of specific beats, Cp/c v 
d im ensionless axial coordinate 
absolute-viscosity coefficient 
kinematic-viscosity coefficient 
dimensionless pressure functions 
density 
Pr/r 

dimensionless compressible temperature function, 
t(z) = (T - Too)/(T w - Tj 

radial shear per unit area on planes parallel to the 
disk 

angular shear per unit area on planes parallel to the 
disk 

dissipation function 
angular disk velocity 

disk surface 


0 disk edge 

o° undisturbed region 


Superscript: 


denotes differentiation with respect to z in eqs. 
(7) and with respect to £ in eqs. (13) 
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APPENDIX B 


COMPARISON OP ENERGY EQUATIONS 

To understand the difference Between the energy equations used in 
reference 4 and herein, the two forms of the energy equation will be 
considered neglecting dissipation for axially symmetric boundary layers. 
These are 

/dT 5r\ 5 ( St\ U 5w\ . . 

pc v M + w 5z) “ 3z 3zJ - p \Sr + r + 3zJ (B1) 

/„ St „ St\ S /.a t\ , („ Sp , „ Sp\ ,_ oS 

pc pV u + w = 3z (^ k 3zj + ( u 3r + w SzJ (B2) 

Equation (Bl) was used in reference 4 and (B2) was used herein. The 
last terms on the right were emitted in both cases so that the remaining 
equations Eire identical except for the different specific heats. Clearly, 
for gases the two equations give different results (as has already been 
pointed out). The question then arises whether the omission of these 
terms is really justified in both cases. If the state equation 

P = fdRT 

is applied to equation (B2), there results 

p (°v - B) ( u H * w H) = 'k ( k H) + BT ( u m + w i§) 

or 

pc v ( u m . + w H) = ^ ( k H) + BT ( u si + w lz) (B3) 

Equation (B3) is identical with (Bl) except that the last terms are in 
a different form. From equation (B3) it is evident that the last term 
is of the same order of magnitude as the convection (left side) term. 
Therefore, with the c v form of the energy equation it is not correct 
to omit the last terms of equation (Bl) for gases as was done in refer- 
ence 4. Looking at (B2), however, it can be seen that within the 
boundary- layer assumptions the last term can be neglected relative to 
the convection terms as was done herein. 



NACA TN 4320 


13 


REFERENCES 

1. von K&rmKn, Th. : Uber laminare und turbulent Reibung . Z.a.M.M. , 

Bd. 1, Heft 4, Aug. 1921,, pp. 233-252. 

2. Cochran, W. G.: The Plow Due to a Rotating Disc. Proc. Cambridge 

Phil. Soc., vol. 30, July 30, 1934, pp. 365-375. 

3. Wagner, C.: Heat Transfer from, a Heated Disk to Ambient Air. Jour. 

Appl. Phys., vol. 19, no. 9, Sept. 1948, pp. 837-839. 

4. Millsaps, Knox, and Pohlhausen, Karl: Heat Transfer by Laminar Flow 

from a Rotating Plate. Jour. Aero. Sci., vol. 19, no. 2, Feb. 1952, 

pp. 120-126. 

5. Ostrach, Simon: An Analysis of Laminar Free-Convection Flow and Heat 

Transfer About a Flat Plate Parallel to the Direction of the Gener- 
ating Body Force . NACA Rep . 111.1, 1953 . (Supersedes NACA TN 2635 . ) 

6. Cobb, E. C., and Saunders, 0. A.: Heat Transfer from a Rotating Disk. 

Proc. Roy. Soc. (London), ser. A, vol. 236, no. 1206, Aug. 2, 1956, 
pp. 343-351. 

7. Young, R. L. : Heat Transfer from a Rotating Plate. 

78, no. 6, Aug. 1956, pp. 1163-1168. 


Trans. ASME, vol 



TABLK I. - BMHSIOHLSBS BOIOTIOJK 


c 

B " 

H " 

a 

G « 

a 

S ' 

3 

t 


— 

H ' 

H 

O ' 

a 

8 ' 

s 

0 

- 1.0205 

0 

0 

- 0.6159 

1.0000 

- 0.3286 

1.0000 

u.o 

0.0001 

- 0.0001 

- 0.8844 

- 0.0001 

- 0,0001 

- 0,0010 

0.0015 

.£ 

-.6676 

-.1673 

-.0179 

-.5987 

.8780 

-.3276 

.9543 

11.2 

.0001 

-.0001 



-.0001 

.0001 

-.0008 

.0013 

.4 

-.5969 

-.2727 

-.0628 

-.5577 

, 7$ai 

-.3285 

.8888 

11.4 

.0001 

-.0001 



-.0000 

.ex 

DOO 

-.0007 

.0012 

.6 

-.2031 

-.3320 

-.1239 

-.5048 

.6557 

-.3222 

.8039 

11.6 

.0001 

-.0001 







-.0005 

.0010 

.© 

-.0635 

-.5578 

-.1954 

-.4478 

.5606 

-.3149 

.7401 

11.8 

.0000 

-.0001 







-.0006 

.0010 

1.0 

.0514 

-.3603 

-.2655 

-.3911 

.4766 

-.3047 

.6781 

12.0 



-.0000 







-.0005 

.0009 

1.2 

.0922 

-.5475 

-.3365 

-.3381 

.4038 

-.2917 

.6184 

12.2 











-.0004 

.0007 

1.4 

.1279 

-.5261 

-.4038 

- .2898 

. 5*11 

-.2768 

.5616 

12.4 











-.0004 

.0006 

1.6 

.1456 

-.2975 

-.4661 

-.2470 

.2875 

-.2598 

.5079 

12.6 





-.6845 





-.0003 

.0005 

1.6 

.1508 

-.2677 

-.5227 

-.2095 

.2119 

-.2419 

.4577 

12.8 











-.0003 

.0005 

2.075 

.1462 

- .2207 

-.5906 

-.1662 

.1905 

-.2166 

.3955 

13.0 





1 






-.0003 

.0004 

2.275 

.1357 

--+1986 

-.6331 

-.1401 

.1599 

-.1984 

.3536 

13.2 











-.0002 

.0004 

2.475 

.1239 

-.1726 

-.6702 

-.1178 

.1342 

-.1807 

.3157 

15.6 











-.0002 

.0003 

2.675 

.1113 

-.1491 

-.7023 

-.0991 

.1125 

-.1639 

.2813 

13.6 











-.0002 

.0003 

2.875 

.0686 

-1881 

-.7300 

-.0832 

.0944 

-.1477 

.2502 

15.8 











-.0002 

.0005 

5,076 

.0665 

-. 109 $ 

-.7537 

-.0698 

.0791 

-.1327 

.2222 

14.0 











-.0001 

.0002 

5.275 

.0752 

-.0934 

-.7740 

-.0585 

.0663 

-.1189 

.1970 

14.2 











-.0001 

.0002 

5 . 47S 

.0650 

-.0794 

-.7912 

-.0491 

.0556 

-.1062 

.1746 

14.4 











-.0001 

.0002 

5.076 

.0659 

-.0674 

-.0009 

-.0412 

.0466 

-.0947 

.1545 

14.6 











-.0001 

.0002 

3.876 

.0478 

-.0570 

-.8183 

-.0343 

.0390 

-.0843 

.1366 

14.8 











-.0001 

.0001 

4.075 

.0408 

-.0482 

-.8288 

- .0269 

.0327 

-.0748 

.1207 

15.0 











-.0001 

.0001 

4 . e75 

.0346 

-.0407 

-.6377 

-.0242 

.0274 

-.0664 

.1066 

15.2 











-.0001 

.0001 

4.473 

.0294 

-.0343 

-.8451 

-.0203 

.0230 

-.0587 

,0939 

15.4 











-.0001 

.0001 

4.675 

.0249 

-.0209 

-.8514 

-.0170 

.0192 

-.0520 

.0626 

15.6 











-.0001 

.0001 

4. 876 

.0210 

-.0243 

-.8567 

-.0143 

.0161 

-.0459 

.0728 

15.8 











-.0000 

.0001 

5.075 

.0177 

-.0204 

-.8612 

-.0118 

.0155 

-.0406 

.0644 

16.0 













.0001 

5.275 

.0149 

-.0172 

-.8649 

-.0100 

.0113 

-.0358 

.0568 

16.2 













.0001 

5.475 

.0126 

-.0144 

-.8681 

-.0084 

.0095 

-.0316 

.0501 

16.4 













.0001 

5 . S75 

.0106 

-.0121 

-.8707 

-.0070 

.0079 

-.0279 

,0441 

16.6 













.0001 

5.875 

5.075 

6.275 

6.475 

.0069 

.0075 

.oofes 

-.0102 

-.8730 

-.8748 

-.8764 

-.8777 

-.0058 

-.0049 

.0067 

-.0246 

.0389 

16.8 

17.0 













.0 

000 

* 0u71 


-.0191 

-.0169 
















:58s 

-.0060 

-.0055 

.0039 

Wtm 

17.4 















6.675 

.0044 

-.0050 

-.8788 

-.0029 

.0033 

-.0149 

,0234 

17.6 















6.875 

.0057 

-.0042 

-.8797 

-.0024 

.0027 

-.0127 

.0206 

17.8 















7.075 

.0031 

-,(XXJ5 

-.8806 

-.0020 

,0023 

-.0115 

.0182 

10.0 















7.275 

.0026 

-.0030 

-.8811 

-.0017 

.0019 

-.0102 

.0180 

18.2 















7.475 

,0028 

-.0025 

-.5617 

-.0014 

.0016 

-.0090 

.0141 

18.4 















7-673 

.0018 

-.0021 

-.8821 

-.0012 

.0013 

-.0079 

.0124 

16.6 















7.575 

.0015 

-.0017 

-.8825 

-.0010 

.001 1 

-.0070 

.0109 

18.6 















8 . 07B 

.0013 

-.0010 

-.Mae 

-.0005 

.0009 

-.0061 

.0096 

10.0 















a. a 

.0012 

-.0013 

-.8530 

-.0006 

.0006 

-.0056 

.0089 

19.2 















a . 4 
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Figure 2. - DimenBionless velocity functions, F => O/kRj G = 7/feE; H = p¥/C^/ 2 (2' , „) 1 / 2 . 
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Figure 3. - Dimensionless temperature distribution. 
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Figure 4. - Ratio of Nusselt number to square root of Reynolds 
number as function of disk angular velocity. 
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